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Abstract: In this work we analyze the effect of smoothing maps containing arrival directions of cosmic rays with a 
gaussian kernel and kernels of the mexican hat wavelets of orders 1, 2 and 3. The analysis is performed by calculating the 
amplification of the signal-to-noise ratio for several background patterns (noise) and different number of events coming 
from a simulated source (signal) for an ideal detector capable of observing the full sky with uniform coverage. We extend 
this analysis for a virtual observatory with two sites, one in the northern hemisphere, the other in the southern, considering 
an acceptance law. 
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1 Introduction 

There are different experiments aiming to study cosmic 
rays H] 12 [3] HJ |5] . They span different energy ranges and 
are located in different positions on the Earth, but they all 
share a common goal: to unveil the origin, propagation and 
mechanisms of acceleration of the cosmic rays, specially at 
the highest energies. 

The identification of possible astrophysical sources and the 
investigation of the magnetic fields which permeates the 
universe are studied by analyzing the arrival directions of 
cosmic rays. The correlation of these directions with the 
large structures in the universe, such as the galactic and su- 
pergalactic planes are considered large scale anisotropies. 
A small scale anisotropy is characterized by the associa- 
tion between arrival directions of cosmic rays and point 
sources, such as stars, distant galaxies and other kind of 
astrophysical objects with angular sizes small enough to 
be considered point-like objects. Since the sources of the 
cosmic rays of the highest energies are still unknown, it is 
an important task to disentangle genuine signals from the 
background. 

2 Wavelets 

Wavelets are localized wave-like oscillating functions be- 
longing to the L 2 space. When operated with a given sig- 
nal, wavelets can be very useful to extract information con- 
cerning this signal, justifying its use for signal processing. 



The continuous wavelet transform (CWT) in two dimen- 
sions may be formally written as: 

$( S) n,r 2 )=/ j f(t,u)*l TltT2 (t)dtdu, (l) 

where s (s > 0, s € R) is the scaling factor and T\ 
and T2 (ji £ R) are the translation parameters. So, the 
CWT decomposes a function f(t, u) in a basis of wavelet 

The function ^s,ti,t 2 (t, u) is obtained by means of scaling 
and translation of a so-called "mother-wavelet" "J: 
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The mexican hat wavelet family (MHWF) and its extension 
on the sphere have been widely used aiming the detection 
of point sources in maps of cosmic microwave background 
(CMB)Q |9], due to the amplification of the signal- 
to-noise ratio when going from the real space to wavelet 
space. 

The MHWF is obtained by successive applications of the 
laplacian operator on the two-dimensional gaussian. A 
generic member of this family is: 
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where tf> is the two-dimensional gaussian {(j){x) - 
J- e ~ x / 2<T ) and the laplacian operator is applied n times. 
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3 Celestial Maps 

Celestial maps are powerful tools to study anisotropics of 
cosmic rays and are ontained by pixelizing the celestial 
sphere taking into account the resolution of the experiment. 
The events map is the celestial map that represents the ar- 
rival directions of cosmic rays in the celestial sphere in a 
given coordinate system. 

Due to limitations of the detector itself, it is impossible 
to determine the exact arrival direction of an event. Each 
event detected is convolved with a probability distribution 
related to the angular resolution of the experiment, this is, 
for each event there is an associated point spreading func- 
tion (PSF). Therefore, it is extremely useful to convolve 
the celestial maps with functions associated to the PSF of 
the detector, aiming to maximize the signal-to-noise ra- 
tio. Mathematically, this process of convolution (filterin^3) 
may be written as 



M f (r Q )=a / M(r)$(r,ro)dfi, 



(4) 



where a is a normalization constant, M(r) is the number of 
cosmic ray events in the direction r, $(r, fa) is the kernel 
of the transformation and fo is the position vector repre- 
senting each point in which the integral is evaluated. In the 
discrete case this process is: 



M f (k) 



(5) 



where M(j) is the number of cosmic rays associated to the 
pixel of index j in the direction fj . 

4 Analysis Procedure 

In the present work we have tested the capability of the 
MHWF to detect point sources embedded in different 
backgrounds. A source was simulated in the position 
(I, &)=(320°, 30°) (galactic coordinates) in the sky, with an 
angular size of 2°. The number of events in the direction 
of the source is 200 and the source is 10% brighter than 
the background. We have also performed a similar analysis 
for a source with the same amplitude and 50 events coming 
from its direction. 

The simulations of the background were performed accord- 
ing to several background patterns: 

• isotropic: isotropic distribution of events; 

• dipole 1: a dipole with excess in the galactic center 
(I, b) = (0°, 0°), with amplitude 7% with respect to 
the background; 

• dipole 2: a dipole with excess in the direction (I, b) 
= (266.5°, -29°), with amplitude 0.5% with respect 
to the background; 



• sources: several sources with different angular 
scales a and amplitudes A in the directions (l,b): 
(0°,0°) [a = 7.0°, A = 100%], (320°, 90°) [a = 
1.5°, A = 5%], (320°, -40°) [a = 0.5°, A = 1%], 
(220°, 10°) [a = 3.0°, A = 5%], (100°, -70°) 
[a = 2°, A = 10%], (240°, 50°) [a = 20°, 
A = 5%], (350°, -80°) [a = 6.0°, A = 0.5%], 
(100°, 50°) [a = 30°, A = 50%], (140°, -40°) 
[a = 4.0°, A = 200%] and (60°, 50°) [a = 3.0°, 
A = 2%]. 

To convolve the filters with the simulated maps we have 
considered two different scenarios for the detector: (1) a 
full sky coverage with a uniform exposition, modulated 
only by the background pattern imposed to the simulations; 
(2) a detector with an acceptance law. In the last case, we 
considered two detectors, one in the northern hemisphere 
with latitude 38° N and longitude 102° W, and the other in 
the souther hemisphere at 36° S and 65° W. The detector 
on the northern hemisphere has an area seven times greater 
than the one in the south, so that the flux of cosmic rays is 
also seven times greater. Also, we have assumed that the 
event are detected according to a zenith angle distribution 
that follows sin9cos9, where 0° < 8 < 60° is the zenith 
angle|[10l. 



In order to verify the power of identification of the 
wavelets, we have calculated the amplification (A) of the 
signal-to-noise ratio, which is given by: 



A = 



(6) 



where Wq is the value of the central pixel associated to the 
source in the non-filtered source map, Wf is the value of the 
same pixel in the filtered source map, <jq is the root mean 
square (RMS) of the non-filtered background map and a f 
is the RMS of the filtered background map. 

According to Gonzalez-Nuevo et al. HI 111 , the maximum 
amplification of the signal-to-noise ratio for a source with 
angular size 70 embedded in a white nois^l background is 
obtained by convolving the maps with a gaussian kernel 
with dispersion 7 = 70- This shall be used as a reference 
to check the consistency of our results and can be seen in 
figure Q] 

The same reasoning can be applied for a source embedded 
in an isotropic background with the acceptance of the de- 
tector. However, as shown in figure [2] the gaussian kernel 
is not the filter that provides the best amplification. 

For each case of signal and background we have plotted 
graphs similar to the ones shown in figures Q] and [2] 

In figure [3] it shown the maximum value of the amplifica- 
tion for each filter for the different background patterns, in 
the case of 800,000 events from the background and 200 



1. In this paper we make no distinction between the processes 
of convolution, smoothing and filtering. We also make no distinc- 
tion between the terms kernel and filter. 

2. White noise has the property of being homogeneous and 
isotropic, with a uniform power spectrum. 
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Figure 1: Amplifications for a source embedded in an 
isotropic background. 
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Figure 2: Amplifications for a source embedded in an 
isotropic background considering the acceptance of the de- 
tector. 

from the source. For all the other cases, the analysis is sim- 
ilar and we can draw the same conclusions. It is interesting 
to notice that the maximum amplification for the MHWF 
kernels are achieved approximately at the same scale, in- 
dependently of the background pattern and the existence of 
an acceptance, whereas the gaussian filter is strongly de- 
pendent on the coverage. 

5 Conclusions 

In this work we have considered a sinQcosQ law for the 
zenith angle distribution for a virtual detector with two 
sites, in both hemispheres, and an ideal detector observ- 
ing the whole sky, with uniform acceptance. We have sim- 
ulated a source with different number of events coming 
from it, embedded in several background patterns. We have 
calculated the amplifications of the signal-to-noise ratio in 
each case and verified that in the case of a source with an- 



gular size 70 embedded in a white noise background, the 
maximum amplification is achieved with a gaussian filter 
with dispersion 7 = 70, as predicted. Also, in all the 
cases with full coverage, the gaussian had a slightly bet- 
ter performance than the MHWF filters. Even though we 
have introduced background patterns which could lead to 
a background with non-white spectral features, such as a 
dipole and several sources, the gaussian kernel provided a 
slightly greater amplification, but only in the case of full 
sky coverage of the detector lfl2ll . 

It is interesting to notice that the existence of an acceptance 
for the detector affects the power of amplification of the 
gaussian filter, independently of the spectral shape of the 
background. The amplification, in this case, is always be- 
low 2, whereas in the case of the MHWF the amplification 
of the signal-to-noise ratio can be greater than 6. 

If the background presents characteristics of non-white 
noise, the maximum amplification of the signal-to-noise 
ratio can be achived by the MHWF kernels, when com- 
pared to the gaussian. Although, when using MHWF, the 
directional information regarding the source is less precise. 
For instance, if we consider a gaussian filter and a white 
noise background, the width of the optimal gaussian filter 
matches exactly with the angular size of the source, thus the 
uncertainty on the position would arise only from the reso- 
lution of the detector rather than the analysis method used. 
By using MHWF, the greater the order of the wavelet, the 
greater will be the difference between the angular sizes of 
the source and the wavelet which provides the maximum 
amplification. Therefore, despite the gain on the power of 
discrimination of the sources, the accuracy on their local- 
ization decreases and this is related only to the method and 
not to the instrument. We conclude that is important to con- 
sider the effects of the acceptance of the detector and their 
impact on the power spectrum of the background in celes- 
tial maps and consequently, on the choice of the best kernel 
to be used. 

The power spectrum of the cosmic rays is not yet fully char- 
acterized, so it is important to have analysis tools that do 
not strongly depend on specific features of the background 
to amplify the signal-to-noise ratio. Since the acceptance 
introduces an unknown spectral shape for the background, 
for analysis involving the whole sky, such as blindsearcl{| 
the MHWF amplifies the signal-to-noise ratio more than 
the gaussian. However, for small scale analysis, smooth- 
ing the maps with a gaussian kernel and reducing the area 
of scan could provide a greater amplification, specially if 
the background within the window of scan can be approxi- 
mated by a white noise with a uniform power spectrum. 

A limitation of this technique is the projection of the ce- 
lestial sphere to the plane. This approximation is good 
enough for small scales, when we use sinO « 0. For 
9 > 7°, however, this approximation introduces bias and 
this method would have to be adapted to work on a spher- 

3. The procedure of blindsearch is the search for astrophysical 
objects within a given window around an observed excess. 
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Figure 3: Maximum amplifications for each filter for all the simulated background patterns. 



ical manifold. Wavelets on the sphere have been used on 
CMB studies lfl3l fT4l [151 [TBI and presented good results. 
The next step of this work is to perform a sistematic study 
relating the effect of the increase in the optimal scale of 
each kernel with respect to the angular size of the source. 
Once this is achieved, a correction factor might be applied 
for a better constraint on the position of the source. Further- 
more, we intend to extend the method to the sphere, which 
allow us to search not only for point sources, but also for 
large scale structures in maps containing arrival directions 
of cosmic rays. 
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